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Abstract. We present a first catalogue of X-ray sources resulting from the central area of the XMM-LSS (Large 
5^ ■ Scale Structure survey). We describe the reduction procedures and the database tools we developed and used to 

| derive a well defined catalogue of X-ray sources. The present catalogue is limited to a sub-sample of 286 sources 

detected at 4<r in the 1 deg 2 area covered by the photometric VVDS (VIRMOS VLT Deep Survey), which allows 
us to provide optical and radio identifications. We also discuss the X-ray properties of a larger X-ray sample of 
536 sources detected at > 4<r in the full 3 deg 2 area of the XMM Medium Deep Survey (XMDS) independently 
of the optical identification. We also derive the logN-logS relationship for a sample of more than one thousand 
I— i ' sources that we discuss in the context of other surveys at similar fluxes. 
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1. Introduction 

The XMDS (XMM Medium Deep Survey) is an X-ray sur- 
vey based on the pooling of a significant fraction of the 
XMM-Newton guaranteed time of three "hardware insti- 
tutes" (IASF Milano for XMM-EPIC, Liege for XMM-OM 
and CEA Saclay for XMM-SSC). 

The XMDS pointings lie at the heart of the full, larger 
XMM Large Scale Structure (XMM-LSS) Survey (see 
Pierre et al., 2004), to which we refer for a discussion of 
the motivations and the choice of the sky field. 

About two thirds of the XMDS area are covered in 
the optical band by the VVDS (VIRMOS VLT Deep 
Survey) both by UBVRI photometry (Le Fevre et al., 
2004a) and by multi-object spectroscopy with VIMOS (Le 

* The full catalogue illustrated in Table 3 is only 
available in electronic form at the CDS via anony- 
mous ftp to cdsarc.u-strasbg.fr (130.79.128.5) or via 
http : // cdsweb . u-strasbg . f r/ cgi-bin/ qcat ? J/A+A/ 



Fevre et al., 2004b), by an associated radio survey at 1.4 
GHz (Bondi et al., 2003) and by a GALEX Early Release 
Observation (Arnouts et al., 2005). 

The XMM-LSS area has been covered by an associated 
radio survey at 74 and 325 MHz (Cohen et al., 2003) and 
will be covered by surveys in other bands like SWIRE 
(Lonsdale et al., 2003), the CFHTLS (see web site 1 ) and 
by our own follow up observations (Pierre et al., 2004). 

Thus, for these areas, multiwavelength information ex- 
ists or will be gathered in the future, allowing investi- 
gations of large scale structures and source populations 
(clusters, galaxies and AGN) at medium redshifts. 

For internal usage within the XMDS and XMM-LSS 
consortia, and in the future also for public access, a 
database has been designed, which is presented in this 
paper together with the first results obtained from the 
analysis of the XMDS data. 



http : //www. cfht .hawaii . edu/Science/CFHLS/ 
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The layout of the paper is as follows. The database tool 
devised to store the results of the XMM-LSS and XMDS 
and associated surveys is described in Section 2. Section 3 
contains the analysis procedure (closely following that of 
Baldi et al., 2002) used to generate source lists from the 
19 XMM pointings of the XMDS using a threshold that 
ensures a small number of spurious sources. Different se- 
lection criteria are then applied to obtain samples used for 
different purposes. A relatively low probability threshold 
is used to compute the logN-logS relation presented and 
discussed in Section 4. The analysis of the X-ray charac- 
teristics is restricted to sources with a signal to noise ra- 
tio larger than 4, and the results are presented in Section 
5. Only a fraction of the XMDS area is covered by the 
VVDS survey. For the latter area the identification pro- 
cedure based on the photometric VVDS data and the de- 
rived optical vs X-ray properties are described in Section 
6. 

Preliminary accounts of this work have been presented 
at recent conferences (Maraschi et al. , 2003; Tajer et al., 
2004). 

2. The database and catalogue 

2.1. The LSS database 

The web site of the XMM-LSS Master Catalogue 2 has 
been designed as a front end to access the database con- 
taining the catalogues, both as a working tool for the 
XMDS and XMM-LSS consortia and, in the near future, 
for public access. 

The query interface uses a Java servlet (communi- 
cating via JDBC protocol to the underlying MySQL 3 
database) to manage the permissions of different groups 
of people to access different (parts of) database tables, to 
perform selections into multiple tables in an user-friendly 
way, to have a quick look at the data and to export results 
as ASCII or FITS files. 

The database presently includes database tables for the 
results of the X-ray pipelines, VVDS photometry for ob- 
jects within a 40" box around an X-ray source, radio and 
other catalogues as reported in Section 6.1; it also pro- 
vides links to data products like X-ray images and expo- 
sure maps, finding charts in the I band, radio maps, etc. 

The database includes also pre-calculated correlation 
tables, which link the sequence identifiers of objects in an 
X-ray table with their neighbours in other tables accord- 
ing to predefined proximity criteria (e.g. a 40" box or a 
6" radius), and greatly speed up queries which involve a 
couple of tables. 

Correlations between more than two tables are man- 
aged similarly by multi-column correlation tables gener- 
ated by our identification procedure described below in 
section 6.1. They are also used to build virtual tables, i.e. 

2 http : //cosmos .mi . iasf . enr . it/~lssadmin/Website/ 
LSS 

3 MySQL is an open source database server developed by 
MySQL AB (see http://www.mysql.com) 



simultaneous views of columns taken from different indi- 
vidual tables, or even of results of algebraic operations 
between different columns. 

Usage of appropriate MySQL syntax during the com- 
position of a query allows the production of results of ar- 
bitrary complexity as output from our database, like e.g. 
ds9 region files, or the DTgX code for the catalogue table 
reported here (Table 3). 

2.2. The catalogue 

In this paper we publish, as a first part of the XMDS cata- 
logue, the results of the identification procedure described 
in section 6.1, i.e. a catalogue of 286 X-ray sources de- 
tected at more than 4ct and located in the VVDS area. 
Further releases will also include sources already present 
in our database and detected at a lower S/N and/or lo- 
cated outside the VVDS area. 

Our current choice allows us to present a catalogue 
inclusive of reliable (though not spectroscopically con- 
firmed) identifications and of the basic optical and radio 
characteristics of the X-ray sources. 

The catalogue published electronically at CDS along 
with the present paper (and available as Online material) 
includes a rather reduced, manageable number of columns. 
A printout sample of the catalogue in tabular form is pre- 
sented in Table 3. Comments on individual sources are 
presented in Appendix A. 

Since our database contains much more information 
than what is included in the published catalogue, we plan 
to open access to our database query engine, in accordance 
with the "legacy policy" of the XMM-LSS survey. 

Usage of our database interface will allow the user 
to make selections on the datasets, to access additional 
columns (including e.g. X-ray fluxes in all bands, photom- 
etry in additional bands, radio fluxes, identification ranks 
and flags, technical information, etc.) which is not prac- 
tical to include here, and to generate (and operate upon) 
expressions involving columns. 

The query interface will be accessible from the time the 
paper is accepted at http://cosmos.mi.iasf.cnr.it/ 
"lssadmin/Website/LSS/Query logging in as user xmds, 
password guest in workspace public. 

3. Observations and data reduction 

3.1. Observations 

The XMDS observations consist of 19 overlapping point- 
ings (typical duration in the range 20-25 ks), covering a 
contiguous area of about 3 deg 2 . The sky location of the 
combined field of view (FOV) of the pointings is shown in 
Fig. 1. 

The observations were performed by XMM-Newton 
(Janscn ct al., 2001), with the EPIC MOS (Turner et al., 
2001) and pn (Struder et al., 2001) cameras, between July 
2001 and January 2003. A basic journal of observations is 
reported in Table 1, containing the ESA dataset identifier 
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Fig. 1. Location of the XMDS pointings on the sky : the ap- 
proximate FOV footprint (with the pn protruding at the actual 
roll angle) is shown (solid lines); the gray colour for field G12 
indicates that it has not been used (see text). The dashed rect- 
angle indicates the area of VVDS photometry. 



that can be used to inquire for further details at the ESA 
XMM archive. 

The XMDS fields are surrounded at present by 32 ad- 
ditional shorter (typically 10 ks) pointings (B01 to B32) 
performed in AOl and A02 as part of the XMM-LSS pro- 
gram. More detailed information is available on an ancil- 
lary page of the XMM-LSS Master Catalogue web site 4 . 

3.2. Generation of the X-ray source list 

The reduction of the XMDS data has been performed in 
Milan using a streamlined version of the pipeline devel- 
oped by Baldi et al. (2002) for the HELLAS2XMM survey. 
We refer to the latter paper for the formalism and proce- 
dure, and only mention here the main points and the few 
differences. 

We used the XMM-Newton Science Analysis System 
(XMM-SAS) v5.4.1, whose tasks emproc and epproc now 
provide a reliable identification of the location of bad pix- 
els. To remove data contaminated by soft proton flares, 
we used fixed thresholds (on the global background above 
10 kcV) of 0.15 counts s" 1 for MOS and 0.35 counts s" 1 
for pn : as a result of this choice, field G12 could not be 
analysed because the background was almost constantly 
above these thresholds. A mosaic of our X-ray images is 
shown in Fig. 2. 

As in Baldi et al. (2002), we ran the source detection 
(inclusive of background computation) and characteriza- 

4 http : //cosmos .mi . iasf . cnr . it/~lssadmin/Website/ 
LSS/Anc/goepic . html 



tion procedure on merged MOS and pn data in five dif- 
ferent energy bands, which we designate as : A (0.3 — 0.5 
keV), B (0.5 - 2 keV), C (2 - 4.5 keV), D (4.5 - 10 keV) 
and CD (2 — 10 keV). We note that sources detected in the 
CD band do not result from the simple sum of detections 
in the separate C and D bands, but from the detection 
algorithm run on the full 2 — 10 keV band, while quanti- 
ties referring to the total band ABCD (0.3 — 10 keV) are 
obtained by summation of the results in the individual 
bands. 

We have used the XMM-SAS tasks eboxdetect and 
emldetect to obtain a list of candidate source positions, 
but net counts and associated parameters are obtained 
from our own program written according to Baldi et al. 
(2002). In particular fluxes were computed for a power law 
spectrum with r=1.7 and 7V#=2.61xl0 20 cur 2 (i.e. the 
average galactic column density in the XMDS field direc- 
tion (Dickey & Lockman, 1990)), using conversion factors 
calculated from response matrices generated consistently 
with the event selections used in our pipeline. In particu- 
lar, as in Baldi et al. (2002), we used a conservative pattern 
selection of single and double events for MOS, and only 
single events for pn. Since the count rates are derived in 
regions including a fixed percentage (68%) of the expected 
flux from a point source (whose radius depends on off-axis 
positions), our fluxes are correct for point-sources. The 
fact that a source may fall close to (or on) an inter-CCD 
gap is taken into account in flux computation by proper 
usage of the EPIC calibration files, and by the usage of 
exposure maps in the measurement of the count rate. This 
should be an acceptable approximation (the same used by 
Baldi et al., 2002). Anyhow, we record in our database a 
flag indicating whether a source is close to a FOV edge or 
to CCD borders or gaps. 

For insertion in our database, we retain only sources 
that have a chance detection probability less than 2x 10 -4 
in at least one energy band. In the 19 useful pointings 
of the survey, we detected 1322 X-ray sources (including 
multiple detections in overlapping regions, see Table 2), 
satisfying the probability threshold P < 2 x 10~ 4 . In this 
work we present subsets chosen with more stringent cri- 
teria: (1) all sources detected with chance probability less 
than 2 x 10~ 5 to derive the logN-logS relationship (sec- 
tion 4); (2) X-ray sources detected with a signal to noise 
ratio larger than 4 in any one of the energy bands defined 
above (where the error erg on the net number of counts is 
calculated using the Poissonian approximation from the 
gross number of counts ctsi mg 

<7 S = 1 + \Jcts img + 0.75 

according to Gehrels (1986) as in Baldi et al. (2002)) to 
study the mean X-ray properties of the brighter end of 
the sample (the 4cr sample, section 5); (3) all sources in 
the 4(7 sample covered by VVDS photometry, for which 
we present the catalogue (section 2.2) and we discuss the 
X-ray/optical properties (section 6). 

The numbers of sources in each sample are reported in 
Table 2. 
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Field 


ESA obs id 


date 


pointing RA 


DEC 


max exposure 2 
(ks) 


shift d ARA 

(") 


ADEC 

(") 


G01 


011268 


0101 


28 Jan 2002 


02:27:25.4 


-04:11:06.4 


24.9 


-1.60 


0.00 


G02 


011268 


0201 


14 Jul 2002 


02:25:54.2 


-04:09:05.6 


9.8 


0.00 


-0.53 


G03 


011268 


0301 


19 Jan 2003 


02:24:45.6 


-04:11:00.8 


21.8 


-0.53 


1.07 


G04 


010952 


0101 


29 Jan 2002 


02:23:25.3 


-04:11:07.6 


25.8 


-1.60 


0.00 


G05 


011268 


0401 


02 Feb 2002 


02:28:05.1 


-04:31:08.1 


23.7 


0.00 


0.00 


G06 


011268 


1301 


26 Jul 2002 


02:26:34.4 


-04:29:00.8 


13.0 


0.00 


0.53 


G07 


011268 


1001 


30 Jan 2002 


U2:2b:2o.o 


C\ A O 1 D7 1 

-04:31:07.1 


23.5 


-0.53 


-0.53 


G08 


011268 


0501 


25 Jul 2002 


ao oo rr /i o 


C\ A OA C\C\ 1 

-04:29:00.1 


18.3 


0.00 


0.00 


G09 


010952 


0601 


01 T OOOO 

31 Jan 2002 


02:22:45.2 


-04:31:11.1 


22.4 


n/a 


n/a 


G10 


010952 


0201 


29 Jan 2002 


02:27:25.4 


-04:51:04.4 


24.8 


-0.54 


0.53 


Gil 


010952 


0301 


02 Feb 2002 


02:26:05.1 


-04:51:06.1 


21.9 


0.00 


0.53 


G12 a 


010952 


0401 


01 Feb 2002 


02:24:45.4 


-04:51:11.2 


n/a 


n/a 


n/a 


G13 


010952 


0501 


03 Jul 2001 


02:23:13.1 


-04:49:03.1 


23.9 


-2.67 


-2.13 


G14 


011268 


0801 


31 Jan 2002 


02:22:04.1 


-04:51:09.7 


13.6 


n/a 


n/a 


G15 


011111 


0101 


06 Jul 2001 


02:27:54.1 


-05:09:02.3 


21.2 


n/a 


n/a 


G16 6 


011111 


0201 


06 Jul 2001 


02:26:34.2 


-05:09:03.1 


3.9 


n/a 


n/a 


G16 b 


011111 


0701 


14 Aug 2001 


02:26:35.2 


-05:08:46.6 


23.7 


n/a 


n/a 


G17 


011111 


0301 


03 Jul 2001 


02:25:14.3 


-05:09:08.4 


22.5 


n/a 


n/a 


G18 


011111 


0401 


03 Jul 2001 


02:23:54.1 


-05:09:09.7 


28.1 


n/a 


n/a 


G19 


011111 


0501 


04 Jul 2001 


02:22:34.0 


-05:09:02.1 


23.7 


n/a 


n/a 



a pointing G12 has not been used because of the high background level (see text) 
b pointing G16 has been repeated since the first instance was curtailed 

c the reported max exposure is the highest value in the exposure maps of the individual cameras, i.e. net of soft proton flares 
d the RA and DEC shifts are the astrometric correction applied (see text) for the fields which overlap the VVDS area and 

for which optical identification has been performed 
e more details available through our ancillary web page http://cosmos.mi.iasf.cnr.it/~lssadmin/Website/LSS/Anc/ 

xmdsepic .html 



Sample 


N tot 


N s 


Ncd 


P < 2 x 10~ 4 


1322 


1166 


419 


P < 2 x 10" 5 


1129 


1028 


328 


4a (detections) 


612 


591 


158 


4cr (independent) 


536 


518 


143 


4a VVDS 


286 


278 


73 



Table 2. Total number of sources and sources detected in the 
B and CD bands respectively for the samples presented in the 
paper. For the P < 2 x 10~ 4 , P < 2 x 10~ 5 and the first 
line of the 4a sample, numbers refer to detections (i.e. include 
multiple detections of the same source in different fields) , while 
for the second line of the 4a sample and for the 4a VVDS 
sample numbers refer to independent sources. 



Note that in the present version we have analysed 
each pointing independently, in spite of the fact that 
the FOVs of adjacent pointings overlap. The same ob- 
ject could therefore be detected in more than one field. 
This fact is recognized only a posteriori (see section 5). In 
what follows we will always specify whether the sample we 
consider includes multiple detections of the same source. 

An independent pipeline is being developed in Saclay 
for the analysis of the entire XMM-LSS, using a wavelet 
technique and therefore best suited for extended sources 



(Pacaud et al. (2005), see a preliminary account in Pierre 
et al. (2004)) and will be used as a basis for the future 
complete XMM-LSS catalogue. 



4. The logN-logS relationship 

We computed the logN-logS distributions in the 0.5 — 2 
and 2 — 10 keV band as follows. 

We have considered all sources with a detection prob- 
ability in the 0.5 — 2 and 2 — 10 keV bands below P = 
2 x 10~ 5 . This is well below our acceptance threshold of 
2 x 10 -4 , ensuring very few spurious sources in the sample 
while giving us a large dataset to compute the logN-logS. 

We have generated a flux limit map for each X-ray 
field and energy band, which contains the faintest flux at 
which a source can be detected at the assumed level of 
significance above the local background. From these we 
have computed the sky coverage plotted in Fig. 3. This 
considers all fields as independent areas, consistently with 
our detection procedure (see section 3.2), so the full sky 
coverage is the sum of all areas surveyed. 

Our faintest flux is ~ 10 -15 erg cm~ 2 s _1 in the 0.5 — 2 
keV band and ~ 7 x 10~ 15 erg cm~ 2 s _1 in the 2 — 10 keV 
band. 
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Fig. 2. Mosaic of our X-ray images in the total (0.3 — 10 keV) band. No exposure map correction has been applied to this figure. 



In Fig. 4 we plot the cumulative logN logS in the 0.5—2 
keV band (upper left panel, 1028 sources) and in the 2 — 10 
keV band (lower left panel, 328 sources). The plots relative 
to individual fields have also been produced and inspected 
and do not show significant deviations from the combined 
distribution. 

The shapes of the (cumulative) logN-logS are clearly 
curved and cannot be represented by a single power law. A 
fit to the cumulative distribution does not allow a mean- 
ingful computation of confidence limits on the parameters, 
because points and errors are not independent (Murdoch 
et al., 1973; Crawford et al., 1970). Therefore we com- 
puted the differential logN-logS by binning (in logarithmic 
space) the number N(F) of sources with flux F into flux 
bins of width AFi (i.e. between fluxes Fj m i n and Fj max ), 
then computing the average sky coverage Ai in the bin. 
The i-th bin of the differential curve n(F) (shown in the 
right panels of Fig. 4) is given by 



N(fj) 
AFiAi 



A statistical error due to the Poissonian error on the num- 
ber of sources is assigned to each bin. Again, we confirm 



that a single power law does not give a good fit. Therefore 
we used a broken power law. 

For the 0.5 — 2 keV band, the best fit is 



n(F) 



6.515 x 10 3 x F- 5 2 - 62 F 15 > 10.58 
384.2 x F^ 142 Fis < 10.58 



All fluxes -F15 are normalized to 10 15 erg cm 2 s 1 . 



90% 



confidence limits on slopes are 2.62_g; 2 2 and 1.421 15 , 
while the break position is in the range 1.06j^Q 22 x 10~ 14 
erg cm~ 2 s _1 . 

For the 2-10 keV band we find 

n(F) = 4.483 x 10 4 x F^ 2 ' 91 

for F15 > 21.4, with 90% errors on slope 2.91±°; 4 ^ and 
break at 2.14±g;|^ x 10~ 14 erg cm~ 2 s" 1 . At lower fluxes 
the slope is not well constrained: the best fit is 1.53, but 
the confidence interval ranges between 0.37 and 2.04. 

In Fig. 4 we show the fit results both in the differential 
and the cumulative forms for a better comparison with the 
literature that gives sometimes cumulative (e.g. Baldi et 
al., 2002; Moretti et al., 2003), sometimes differential plots 
(e.g. Yang et al., 2004; Harrison et al., 2003). Given the 
uncertainties in the fits, we prefer the comparison with the 
actual data where possible. 
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Fig. 4. LogN-logS of sources detected at P < 2 x 1CP 5 in band B in the cumulative (upper left panel) and differential (upper 
right panel) form and the same in band CD (lower left and right panels). The units of n(F) are number per 10 -15 erg cm -2 s _1 
and deg -2 . Our best fits are reported as solid lines on the differential plots and have also been converted to the cumulative 
form. 



Our probability threshold P < 2 x 10~ 5 is one of 
those used in Baldi et al. (2002). Since our procedure 
is very similar to theirs and the flux range is nearly the 
same, the logN-logS should in principle be almost iden- 
tical. We find consistency (within errors) with the inte- 
gral HELLAS2XMM distributions for the lower fluxes: 
fo.5-2 £ 3 x 10~ 14 erg cm" 2 s" 1 and F 2 - W <, 2 x 10~ 14 
erg cm~ 2 s -1 , while at higher fluxes we are nominally 
lower, but, taking into account the purely statistical error 
on the number of sources, we are consistent with the lower 
envelope of the Baldi et al. (2002) relationship. 



A good agreement is also found with the Chandra 
Large Synoptic X-ray Survey (CLASXS, Yang ct al., 2004) 
and with the Serendipitous Extragalactic X-ray Source 
Identification (SEXSI, Harrison et al., 2003) program, in 
the common flux range. We compared the differential plots 
and found consistency with the CLASXS survey in the 
soft band (no plot in that band is given in Harrison et 
al., 2003). In the hard band we are consistent with the 
SEXSI survey, but we are systematically lower than the 
CLASXS points; this is not surprising, given that Yang et 
al. (2004) find that their total counts at F 2 _ 8 ~ 10~ 14 erg 
cm~ 2 s _1 arc ~ 70% higher than those of Moretti et al. 
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Fig. 3. Sky coverage for sources detected at P < 2 x 10 -5 
in the B band (solid line) and in the CD band (dashed line) 
computed considering all fields as independent (see text). 

(2003), who use Baldi et al. (2002) data for the flux range 
we are interested in. 

The sky location of the XMDS survey has been in- 
tentionally chosen in order to avoid known bright X-ray 
sources and this could therefore explain, at least partially, 
the low density observed in our survey for fluxes above 
a few 10~ 14 erg cm" 2 s _1 . On the other hand, Yang 
et al. (2004) interpret their overdensity as an indication 
of an underlying large scale structure. We note that the 
CLASXS, like the XMDS, is constructed on contiguous 
pointings, while the HELLAS2XMM and the SEXSI sur- 
veys are serendipitous surveys: it seems reasonable that 
the number count difference we find between the different 
surveys is caused by cosmic variance. 

5. X-ray properties of the total XMDS 4<r sample 

One of our criteria for the definition of a sample is the 
presence of a localized X-ray count excess with signal to 
noise ratio (as defined in section 3.2) > 4 in at least one 
of the energy bands (the "4<r sample"). This results in 
612 detections. Taking into account multiple detections 
(same source detected in different overlapping fields) the 
total number of independent sources reduces to 536. All de- 
tections in different pointings with X-ray positions closer 
than 6" have been automatically considered as multiple 
detections of the same source. For larger separations (up 
to 18", which can be considered as a boundary for sus- 
picious overlaps, since there are only 2 pairs of sources 
detected in the same field which are closer than such a 
distance, one at 8.7" and one at 15.9"), we examined each 
of the 8 pairs and recovered 2 cases: one is included in 
our catalogue and is discussed in Appendix A, the other 



has a distance of 6.4", but one detection is outside the 
FOV of one camera. For all multiple detections we used 
the position and derived parameters from the detection 
that gave the highest significance, unless the source falls 
outside the FOV of one camera, or close to an inter-CCD 
gap. We do not detect any significant flux variation in the 
multiply detected sources with the possible exception of 
three, which are however close to inter-CCD gaps. The 
apparent flux decrease in the 0.5 — 2 keV band from the 
brightest observation in each couple would be of -35%, 
-65% and -28%, with a significance of ~ 3.5, 3.9 and 4.5 
a respectively. 

5.1. Fluxes and hardness ratios 

The flux distribution of the sources in the 4ct sample is 
shown in Fig. 5 for sources detected in the energy bands 
0.5—2 (left panel) and 2—10 keV (right panel) respectively. 

We find that 518 sources (~ 97%) are included in the 
sample because of their detection in the 0.5 — 2 keV (B) 
band, while only 143 (~ 27%) satisfy the detection cri- 
terion in the 2—10 keV (CD) band. Clearly most of the 
latter are also detected in the B band. We do not have 
any source included in the sample because of its detection 
in the 4.5 — 10 keV band only, and only 16 hard sources 
(~ 10%) arc not detected in the 0.5 — 2 keV band. As 
indicated by Fig. 5, we find a posteriori Fx ~ 2.5 x 10~ 15 
erg cm" 2 s _1 as the flux limit in the B band and Fx 
~ 2 x 10~ 14 erg cm" 2 s _1 in the CD band. 

A spectral analysis of the brighter sources is in 
progress. Here we present the results of the hardness ratio 
analysis that can be done for a significantly larger num- 
ber of sources, and is the only spectral analysis possible 
for the fainter sources in the sample. 

To this end we note that, even if a source does not 
satisfy the detection threshold in a band different from 
the one that motivated the inclusion in the sample, we 
have an estimate of the counts and fluxes in all bands 
considered, which we can use to get spectral information. 
Let us define: 



HR 



CR 



2-10 



CRq. 



5-2 



CR2-IO + Ci?o.5-2 



HR c b — 



CR 



2-4.5 



CRq, 



5-2 



Ci?2-4.5 + CRo.5-2 



CR. 



4.5-10 



CR 



2-4.5 



Ci?4.5_io + CR2-4. 1 , 



where CR are the net count rates in the given energy band 
corrected for PSF and vignetting. We only include sources 
for which we have a reasonable estimate of the flux (i.e. a 
signal-to-noise ratio > 2) in the 4.5 — 10 keV band : upper 
limits in such a band do not carry significant information 
since they provide only an upper limit to HR<[ C , which is 
expected for fainter sources with the canonical unabsorbed 
power law spectrum. On the contrary upper limits from 
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Fig. 5. Flux distribution in band B (left panel) and in band CD (right panel) for the total 4a sample. The hatched area gives 
the distribution of the sources detected at > 4a in the concerned band. 
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Fig. 6. Hardness ratio HR vs 2 — 10 keV flux for all sources 
detected at > 4a in the CD band (asterisks). Only a few typical 
error bars are shown to avoid figure overcrowding, inclusive of 
all those detected at < 2a in the B band. Dashed lines indicate 
hardness ratios expected for an absorbed power law spectrum 
with energy index F = 1.7 and intrinsic Nh as reported in the 
label. 



the 0.5 — 2 keV band are interesting as they provide lower 
limits to HRci,. 

We verified that the values obtained from the com- 
bined count rates from different cameras are consistent 
with those computed from single camera data. 



In Fig. 6 we plot the main hardness ratio HR as a 
function of the 2 — 10 keV flux. Also shown in the figure 
are dashed lines indicating the hardness ratios expected 
for a simple power law model with spectral index T = 1.7 
and increasing values of the absorbing column Nh- This 
plot shows that ~ 30% of sources detected in the CD band 
have hardness ratios consistent with Nh > 10 215 cm~ 2 , 
for z = 0. Note that if sources are at higher redshifts, 
the horizontal lines for different values of Nh shift down- 
wards, i.e. the same hardness ratio corresponds to a higher 
absorption column. 

Within our current sample, we do not hnd a spectral 
hardening with decreasing X-ray flux as reported by other 
authors; however our statistics in the hard band is lower 
than those from other surveys. On the other hand, it ap- 
pears from numbers in Table 2 that the fraction of sources 
in the hard band increases when the sample is larger (from 
26% in the 4cr sample to 32% in the total P < 2 x 10~ 4 
sample). Since a less restrictive threshold implies a lower 
flux limit, we can infer that the fraction of hard sources 
should increase as the flux decreases. 

Fig. 7 shows the so-called X-ray colour-colour diagram, 
i.e. the plot of the hardness ratio between the harder bands 
HRdc vs the hardness ratio between soft and hard bands 
HR c b- 

Sources with HR consistent with Nh > 10 22 cm~ 2 in 
Fig. 6 occupy the same locus in Fig. 7 (or Nh > 10 23 cm~ 2 
if we consider that sources are at z ~ 1). This singles them 
out as promising candidates for more detailed studies. 

In spite of the relatively small number of objects, it is 
worth noting that a few sources with low values of HR c b 
and high values of HRd c are present: these are not consis- 
tent with a simple power law model. In particular, their 
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Fig. 7. Hardness ratio HRdc vs HR c b for sources in the > 4u sample (left panel), which are also above 2a in band D. Typical 
error bars are shown (dashed when the significance in the involved bands is below 4a). Grids in right panel represent hardness 
ratios expected with a simple power law model for values of spectral index and Nh reported in the labels at z = (empty 
circles and dashed lines) and at z = 1 (full squares and dotted lines) . 



location in the diagram at HRdc > —0.2 and HR c b < —0.5 
could suggest a concave spectral shape, due to the rela- 
tive lack of signal around 2 keV. Such a shape could result 
from the superposition of an unabsorbed component on 
a flat/absorbed power law as often observed in obscured 
AGNs, where the unabsorbed component can be due to 
scattered or transmitted AGN continuum or to thermal 
and non thermal emission associated with starburst re- 
gions. 

6. The XMDS/VVDS 4a sample 

6.1. The identification procedure 

For the generation of our sample of candidate identified 
X-ray sources we have used the following data sources 
(database tables and data products). 

- The database table of X-ray sources from our pipeline. 

- The optical tables with the VVDS UBVRI(JK) pho- 
tometry including both sources with reliable photo- 
metric data and objects whose photometry is flagged as 
unreliable because the objects are saturated, masked or 
close to the borders of the VVDS area. In addition, we 
have used the associated finding charts. Data quality 
and limiting magnitudes are discussed in McCracken 
ct al. (2003) and Radovich et al. (2004). 

- The table with the VIRMOS 1.4GHz radio catalogue 
(Bondi ct al., 2003) which covers the 1 deg 2 VVDS 
area at a depth of 80 /xJy and a resolution of 6". 

- The table with the XMM-LSS radio catalogue (Cohen 
et al., 2003), a shallow survey at 74 (4 mJy, 5.6 deg 2 ) 
and 325 MHz (275 mJy, 110 deg 2 ) over a larger area 



encompassing the XMM-LSS fields. Occasionally we 
used the associated radio maps. 

— The table with the tentative cluster candidates 

— The NED and SIMBAD database tables, containing 
pointers to objects within 20" from an X-ray source. 
It should be noted that the VIRMOS 1.4GHz cat- 
alogue has recently been inserted also in SIMBAD, 
which could result in a redundant association. 

The identification process has followed a multi-step 
procedure, that also involved the determination of reli- 
able X-ray/optical angular separation criteria for identifi- 
cations and the testing of more automated identification 
procedures for future works. 

First of all, we used our database to create the sam- 
ple of X-ray sources that are within the VVDS area (i.e. 
for which we can provide a FITS finding chart, covering 
a 40" box around the X-ray position, that we visually in- 
spected). 

For each unique X-ray source we have then conser- 
vatively considered all optical objects within a radius of 
6" from the X-ray position (see e.g. Waskett et al., 2004). 
This radius is considerably larger than the nominal SAS 
(emldetect) position error combined with the 2" XMM 
astrometrical uncertainty. For such objects we have posi- 
tions, magnitudes and a crude optical classification (star- 
like vs extended) from the VVDS database. 

We have identified the optical counterparts according 
to the following criteria: a) unique optical counterpart [67 
objects]; b) closest and brightest (Amag>2) counterpart 
[166 objects]; c) no optical counterpart down to I~25, the 
magnitude limit of McCracken et al. (2003) [5 objects, of 
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Fig. 8. Examples from our identification procedure : some representative I band finding charts (size 40"x40") with overlays of 
the X-ray error circles (white thick at 4" radius), the positions of all optical sources within 6" (cyan X's and crosses respectively 
for point-like and possibly extended objects), the positions of NED catalogue objects (yellow diamonds) and the radio source 
error box (green rectangles) or extent (green ellipses). The printed figure is in grayscale : see the online edition for a colour 
version. Three cases are shown : an unique counterpart, an identification preferring a radio source, an ambiguous identification. 



which 3 with no counterpart at all, 2 with a radio coun- 
terpart]. We also used the presence of a radio counter- 
part as a preference criterion. The remaining cases [55] 
do not allow us to assign an unambiguous optical coun- 
terpart until spectroscopic confirmation of the source is 
available, so all possible optical counterparts are retained. 
22 cases were flagged as suspect groups, clusters or simply 
crowded fields, of which half have however a dominant, 
unique counterpart. Examples of finding charts used in 
our identification procedure are shown in Fig. 8. 

After a first pass with this procedure, we found evi- 
dence in some X-ray fields of a systematic rigid shift be- 
tween the position of the X-ray sources and their coun- 
terparts (as high as 3" in field G13). We profitably used 
the pairs identified at such stage to compute (via the SAS 
task eposcorr) a rigid shift correction (see Table 1), and 
inserted in the database two additional columns with cor- 
rected X-ray positions. No rotational correction was re- 
quired. See also section 6.2 and Fig. 9. After this astro- 
metric correction we repeated our inspection and ranking 
procedure. 

This time-consuming identification process involved a 
number of iterations and intermediate steps, which we 
have recorded in the database with different flags and 
fields (that could identify for example the presence of a 
visually crowded field, the cross-correlation with a detec- 
tion of the same source in a different field, the source and 
quality of the optical photometry). These will remain as 
part of the record for each source for future automatic 
database queries, but are not included in the catalogue 
attached to this paper. 

The final identifications are reported in the catalogue 
presented here and have associated numerical codes that 
describe the quality (rank) of the identification, as de- 
scribed in point M of Table 3. 

As a final step, we have computed the probability of 
chance coincidence between the X-ray source and the op- 
tical candidates within the given radius according to the 



formula by Downes et al. (1986) used e.g. by Brusa et al. 
(2003) or Mignoli et al. (2004) 

p = 1 — exp(—n n(< m) r 2 ) 

where r is the distance between the X-ray source and the 
concerned optical candidate, and the density n(< m) of 
optical objects having I magnitude brighter than the mag- 
nitude m of the candidate counterpart was computed us- 
ing the full VVDS catalogue (McCracken et al., 2003). A 
p < 0.01 criterion (similar to the one used by Brusa et 
al. (2003) and more conservative than the one used by 
Mignoli et al. (2004), corresponding in our case e.g. to the 
probability of finding an 1=21 object at 2") applied to the 
X-ray/optical pairs produces a list of identifications very 
similar to our visually identified sources. Namely 93% of 
our best identifications (rank 0) and 84% of the rank 1 
associations have p < 0.01 (and only one of the rank 
identification has p > 0.02), while only 21% of the am- 
biguous (rank 2) cases have p < 0.01. 

At the present time, spectroscopy is in progress within 
the framework of different programs, and only a few of our 
identifications are spectroscopically confirmed. 

6.2. Optical vs X-ray properties 

Fig. 9 shows the angular separation between the position 
of the best optical counterpart identified according to the 
positional and photometric criteria (as described in section 
6.1) and the position of the X-ray source (after astrometric 
correction). 

One can note the absence of systematic shifts. We also 
verified that there is no correlation between the X-ray to 
optical distance and the off-axis position in the field of 
view, while there is a mild correlation with count rate, in 
the sense that sources with a large X-ray to optical sepa- 
ration are among the weakest. We note moreover that 93% 
of the optical counterparts are within a 2" radius from the 
X-ray source (99% within 4"). We therefore conclude that, 
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Fig. 9. Distance in RA and DEC between the X-ray source 
position (inclusive of astrometric correction) and the optical 
counterpart for well identified sources (ranks and 1 in the 
database). Two fiducial radii of 2" and 4" are shown. 
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Fig. 10. R magnitude distribution for well identified sources 
(ranks and 1 in the database, solid histogram) and for optical 
objects within 6" from the ambiguously identified sources (rank 
2, dashed histogram) 



after astrometric corrections have been applied, the search 
radius for optical counterparts can be safely reduced to 4". 

Figure 10 compares the distribution of optical mag- 
nitudes for the identified sources and for the ambiguous 
identifications, i.e. X-ray sources for which there are 2 
or more candidate identifications in the search radius. 
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Fig. 11. X-ray to optical flux ratio (band CD vs R magnitude) 
as a function of the X-ray flux (in band CD). The sources shown 
are those in the 4a VVDS sample, detected above 2a in band 
CD. Upper limits (on the X-ray to optical ratio) correspond 
to sources photometrically saturated (optical magnitude likely 
underestimated). Lower limits corresponds to sources without 
optical counterpart (1 association with a radio source only). 
Open symbols correspond to rank 2 sources (ambiguous iden- 
tification, more than a counterpart, of which the brightest is 
used to derive the optical flux). Diagonal lines indicate loci of 
constant R magnitude while horizontal dotted lines mark the 
region of canonical AGNs (0.1 < F X /F R < 10) 



As expected, unidentified sources populate the region of 
fainter magnitudes where the density of possible optical 
counterparts makes chance coincidences more probable. 
Spectroscopy or at least multicolor photometry is needed 
to increase the number of identifications and to deter- 
mine the redshifts and types of the candidate counter- 
parts. Work on the optical and multifrequency informa- 
tion already acquired within the VIMOS and XMM-LSS 
consortia is in progress. 

We compute the X-ray to optical flux ratio using our R 
magnitudes (in the AB system) according to the formula 

log{F x /F R ) = logF x + 5.51 + Rab/2.5 

where the constant derives from the standard conversion 
of AB magnitudes into monochromatic flux and from the 
integration of the monochromatic flux over the filter band- 
width assuming a flat spectrum. 

The X-ray to optical ratio for our identified sources is 
shown in Fig. 11 as a function of the X-ray flux (2 — 10 
keV band). We have also included the ambiguous identi- 
fications (marked with different symbols) assuming that 
the brighter object is the true counterpart. The latter val- 
ues should therefore be considered as lower limits for the 
X-ray to optical flux ratio: a few objects that are now 
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Fig. 12. X-ray (2 — 10 keV) to optical flux ratio vs hardness 
ratio for sources of the 4<r VVDS sample. Only sources with a 
significance of at least 2a in one band (and of course 4<r in the 
other) are shown, provided they have a rank 0-1 identification 
(one case has only a radio counterpart). Upper limits are due 
to "saturated" optical counterparts which only yield a lower 
limit to the optical flux. 



in the "normal AGN" region could migrate upwards and 
populate the obscured AGN region. 

As expected, most (~ 80%) of the sources lie in the 
region of classical AGNs (0.1 < F X /F R < 10), - 3% are 
in the starburst - normal galaxy range (Fx/Fr < 0.1) 
and ~ 17% (including optically blank fields) in the region 
expected for optically obscured AGNs (Fx/Fr > 10). 

The percentage of high X-ray to optical ratio sources is 
only slightly lower than in other X-ray surveys of similar 
depth (e.g. 20% in Fiore et al., 2003). However, since the 
4cr VVDS sample is not fully identified and since our iden- 
tification procedure favours optically brighter objects, the 
real fraction of sources with Fx/Fr > 10 could be higher. 

Fig. 12 shows the X-ray to optical ratio vs the hard- 
ness ratio. As shown recently by Delia Ceca et al. (2004) 
this diagram is a useful tool to separate stars from the 
extragalactic population, and to locate different kinds of 
AGN. Since we have used the X-ray flux in the 2 — 10 
keV band (for comparison with Fig. 11), stars are not in- 
cluded in our diagram as they are not detected in the 2 — 10 
keV band. We verified however that our confirmed stars 
(via SIMBAD or via preliminary results of VIMOS spec- 
troscopy) fall in the expected location in the ^0.5-4.5/^1 
vs HRcb plot. Most of our sources populate the expected 
locus of AGNs, with a number of objects with very high 
(10-100) X-ray to optical ratios (suggesting optical extinc- 
tion) and a few with hard X-ray spectra, pointing to X- 
ray absorbed nuclei. In a forthcoming paper we will take 
advantage of the full multifrcquency and spectroscopic in- 



formation to discuss the properties of the populations se- 
lected by the X-ray survey with respect to those of the 
optical survey. 
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Table 3. Sample of our X-ray catalogue. We report here only the first few lines of the catalogue to show the layout. The full 
catalogue will be available at CDS, as "Online material" or through our database as described in Section 2.2. Comments on 
peculiar sources are reported in Appendix A. 
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A internal sequence identifier 

B IAU catalogue name of the object (built from uncorrected coordinates as provided by the XMM-SAS) ; the X-ray 
properties of a source, i.e. columns from B to L, appear once only on the first catalogue entry for such a source. 
If the X-ray source has further possible counterparts, these columns are blank and their value is the same as the 
first entry with same sequence identifier. 

C astrometrically corrected RA (J2000) 

D astrometrically corrected Declination (J2000) 

E number of net counts in band B 

F error on net counts in band B 

G number of net counts in band CD 

H error on net counts in band CD 
I flux in 1CP 15 erg cm~ 2 s _1 in band B 

J flux in 1CP 15 erg cm~ 2 s _1 in band CD 

K band in which S/N ratio is the highest 

L gap flag : can assume a value of (good), 1 (fair), or 2 (poor) related to the fact that the X-ray source falls on 
inter-CCD gaps in one or more of the EPIC cameras, or close to the edge or outside the FOV of one EPIC camera 
and is based on the ratio of the punctual exposures in each camera to the maximum exposure in the FOV also for 
each camera. It can also assume the value -1 if the X-ray source is outside the FOV of at least one camera (there 
are no such cases in the 4ct sample). 

M identification rank : gives the quality of the identification with the counterpart, coded as optimal e.g. only one 
counterpart, or brightest and closest; 1 good; 2 (dubious or ambiguous, e.g. typically more than one possible 
counterpart listed); 4 unidentified source 

N distance in arcsec between X-ray source and counterpart in the following order of preference : the VVDS optical 
counterpart if any; or the radio counterpart in the VIRMOS 1.4GHz catalogue (Bondi et al., 2003); or the radio 
counterpart in the XMM-LSS 325 MHz catalogue (Cohen et al., 2003); or the catalogued NED object; or the 
catalogued SIMBAD object. It is blank if there are no counterparts. 

O R magnitude of counterpart from the VVDS catalogue. R and I magnitudes are non-blank when there is a valid 
measurement in the VVDS catalogue. They may be left blank also in case of radio counterparts only, or of (usually 
very bright) objects present only in external (SIMBAD) catalogues. 

P error on R magnitude 

Q I magnitude of counterpart 

R error on I magnitude 

S photometry reliability flag (derived from McCracken et al., 2003) : good; 1 flagged as possibly bad, but not 

saturated; 2 flagged as bad because saturated 
T a V indicates presence of a radio counterpart in VIRMOS 1.4GHz catalogue (Bondi et al., 2003) 
U an L indicates presence of a radio counterpart in the XMM-LSS 325 MHz catalogue (Cohen et al., 2003) 
V name of counterpart, precedence order is NED name, then SIMBAD name including VIRMOS 1.4GHz catalogue 

then name in the radio catalogue by Cohen et al. (2003) 
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Appendix A: Comments on individual sources 

- XMDS J022315.4-042558: it is identified with a bright 
spiral, present also in the NED database (b_J=19.75), 
which however falls in a bad photometric area (close 
to VVDS border) and has no valid optical magnitudes. 

- XMDS J022316.3-042007: it has as counterpart the 
bright (V=11.0) star BD-04 388 which also falls in 
a bad photometric area close to VVDS border, and 
has no valid optical magnitudes ; even if it didn't, it 
would have been saturated. The distance reported in 
the catalogue is overestimated, since the bright source 
in the VVDS catalogue is split in two by a saturated 
spike. 

- XMDS J022330.9-044235: it is considered a dubious 
(rank 2) identification, since the closest object reported 
as rank 2 is faint (1=24.5) and at 5.4". It is a borderline 
case, which could perhaps be a blank field. 

- XMDS J022337.3-040938: it coincides with a bright 
saturated (I < 14.6) starlike object, which matches a 
7.9 mJy 325 MHz radio source, which is a lobe com- 
ponent of the multiple source CRK2003 J0223.6-0409. 
Another component is instead close to the NED object 
NVSS J022336-040932. Could it be that the star is a 
foreground object with respect to the radio source ? 
The X-ray optical distance reported in the catalogue 
is overestimated, since the bright source in the VVDS 
catalogue is split in two by a saturated spike. 

- XMDS J022403. 7-043303: it has as counterpart the 
16th magnitude object Markarian 1036, which is un- 
fortunately saturated in the VVDS. 

- XMDS J022403.9-041326: it falls inside a group of 
weak uncatalogued objects, which are affected by the 
spike of a bright star situated north of them. A ra- 
dio object detected both in the VIRMOS 1.4GHz 
catalogue and in our own observations at 325 MHz 
(CRK2003 J0224.0-0413) is at 13.6" south. This 
source has a duplicate detection in a different pointing 
at 11.7", and both are inside an unpublished galaxy 
cluster (Andreon, private communication). 

- XMDS J022413. 2-045723: it falls at a blank position 
in a small group of objects, unconfirmed yet as a spec- 
troscopic group, of which the closest is at 2.6". Note 
that this source is in field G18, which is not yet astro- 
metrically corrected. 

- XMDS J022445.4-043656: it is considered unidentified 
(rank 4), falling in a blank field since the closest object 
(rejected) is faint (1=25.7) and at 5.0". 

- XMDS J022450.9-042903: its brightest and closest 
counterpart is a 21 mag object at 0.5", however 
at 2.8" there is a radio source detected in the 
VIRMOS 1.4GHz survey as well as in our own 325 
MHz data (CRK2003 J0224.8-0429).We kept it as a 
rank 2 source. 

- XMDS J022454.7-040628: its brightest and closest 
counterpart is a 21 mag object at 1.1", however there 
is a fainter object at 2" which has also 1.4 GHz radio 
emission. 



- XMDS J022456.0-041725: it appears to match posi- 
tionally an uncatalogued enhancement in the periph- 
ery of the NED spiral galaxy 2MASX J02245534- 
0417351 (distance of the nucleus from the X-ray source 
14.8"). 

- XMDS J022524.8-044042: it is associated with a satu- 
rated bright elliptical, also listed in NED (b_J=20.04) 
and as a VIRMOS 1.4GHz radio source. 

- XMDS J022530.6-041418 is one of the sources men- 
tioned in section 5 for which the observations in 
two different pointings give a marginal flux difference 
(3.5a). Its counterpart is a saturated object (nominal 
VVDS 1=16.21), corresponding to the NED elliptical 
galaxy 2MASX J02253077-0414183. 

- XMDS J022544.6-041936 is the other source in the 
VVDS/4er catalogue mentioned in section 5 for which 
the observations in two different pointings give a 
marginal flux difference (3.9ct). Its counterpart is also 
a saturated object (nominal VVDS 1=17.49), corre- 
sponding to the NED elliptical galaxy APMUKS(BJ) 
B022313.84-043304.3. 

- XMDS J022606.6-040259: it is assigned rank 2 because 
the likely counterpart is not too close to the X-ray po- 
sition and appears to be an unresolved double object. 

- XMDS J022606.9-043316: its counterpart should be 
a background object masked behind the bright 
(V=10.37) star BD-05 456 (which is too far to be the 
counterpart). For this reason, the line-of-sight associ- 
ation is classified as rank 2. 

- XMDS J022713. 6-043910: it has as counterpart the 
bright (V=10.12) star BD-05 461 not included in the 
VVDS catalogue; if it were, it would have been satu- 
rated. 

- XMDS J022714.2-042645: the finding chart is satu- 
rated by a bright star, BD-05 480, located southwards 
with respect to our target, whose position is partially 
behind a spike from the star. For this reason there are 
no objects listed in the VVDS catalogue, although an 
object is clearly visible behind the spike at approxi- 
mately RA=02:27:14.16 DEC=-4:26:45.0. 

- XMDS J022721. 7-044153: it has no nearby optical 
counterparts, but a very good positional coincidence 
(0.8") with a VIRMOS 1.4GHz radio source. Note that 
a 17 mag saturated object is however present 6" north. 

- XMDS J022725. 1-043656: it is considered unidentified 
(rank 4), falling in a blank field since the closest object 
(rejected) is faint (1=25.8) and at 4.2". 

- XMDS J022735.6-041317: it corresponds to a plethora 
of very faint objects all located inside an extended ra- 
dio halo, whose central position matches very well the 
X-ray position. 

- XMDS J022735. 7-041122: there are two possible opti- 
cal counterparts which are both located at the periph- 
ery, but inside the nominal extent, of the extended 
325 MHz radio source CRK2003 J0227.6-0411 cen- 
tered 19" to the north-east. A 1.4 GHz point-like 
source,which is a component of a terribly elongated 
multiple object, is instead located in the X-ray error 
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radius between the two optical objects, while the NED 
object NVSS J022736-041117 is located at 9" mid- 
way between the X-ray source and the 325 MHz radio 
source. 

- XMDS J022742. 1-043607: there are no optical cata- 
logued sources; an unresolved double source is present 
at 4" just outside of the nominal error radius of the 
X-ray source and within 1.1" there is a 1.4 GHz radio 
point source. The finding chart is affected by the spikes 
of a nearby bright star. 

- XMDS J022748.2-041013: it has no optical counter- 
part at all, hence is flagged as rank 4, blank field. 

- XMDS J022812. 5-045717: it is in field G15 which has 
not yet been astrometrically corrected. 

Concerning spectral peculiarities, we can use as an indi- 
cator the fact that the energy band with the best signal 
to noise (column K of the catalogue) is not band B. We 
have only one source ( XMDS J022410. 7-042759) which is 
very soft (detected only in band A), and 10 cases with the 
best signal to noise in band CD (of these 7 are among the 
16 not detected in the softer bands mentioned in Section 
5.1). 
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We report here the complete file with our catalogue. For an explanation of the columns refer to the sample table 
3 in the main text. 
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